Fidelity enhancement by logical qubit encoding 
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Abstract 

' We demonstrate coherent control of two logical qubits encoded in a decoherence free subspace 



(DFS) of four dipolar-coupled protons in an NMR quantum information processor. A pseudo-pure 
fiducial state is created in the DFS, and a unitary logical qubit entangling operator evolves the 
system to a logical Bell state. The four-spin molecule is partially aligned by a liquid crystal solvent, 



CN ■ which introduces strong dipolar couplings among the spins. Although the system Hamiltonian is 

> ■ 

never fully specified, we demonstrate high fidelity control over the logical degrees of freedom. In 
fact, the DFS encoding leads to higher fidelity control than is available in the full four-spin Hilbert 
' space. 
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The successful transition from quantum information processing to quantum computation 
will require the ability to efficiently control qubits in the presence of noise. Decoherence 
free subspaces (DFSs) are some of the most efficient schemes of avoiding decoherence from 

n p n 

noise sources with underlying symmetries [1|, 12, 13[. There have been many experimental 
demonstrations of coherent control in DFSs including demonstrations of multiple qubit con- 
t.o, flflg.flay.yaQQ. noweve, . e^a^pl. thee„cod.n. „ot protect 
against the physically relevant noise sources, and the available control elements did not di- 
rectly match the generators of logical control. Here we explore a DPS encoding of strongly 
dipolar-coupled spins where the logical encoding quite naturally fits the internal Hamilto- 
nian structure. This work builds on extensive theoretical investigations of DPS encoding 



16|, 



17(1 . To achieve a similar 



for systems with time-dependent exchange couplings Is], 
internal Hamiltonian structure we employ liquid crystal solvents to partially align the spin 
system and to reintroduce the spin-spin dipolar coupling jl8| . 

Liquid crystals in a strong external magnetic field are partially ordered. This partial 
ordering restricts the thermal motion of molecules dissolved in the liquid crystal material, 
and consequently, the solute molecules have a preferred orientation, and the orientationally- 
dependent intramolecular dipolar interactions do not average to zero [18]. However the 
translational motion of solute molecules is not restricted, and intermolecular dipolar cou- 
plings do average to zero. The resulting system Hamiltonian for a liquid crystal solvent 
system is 
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(1) 
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where Vj is the resonance frequency of the j*^ spin, djk is the dipolar coupling strength 
between spins j and k, Jjk is the corresponding scalar coupling strength, and the sums are 
restricted to spins within the molecule. In liquid crystal solvent systems, dipolar coupling 
strengths can reach several kHz, and for a given pair of spins the dipolar coupling is typically 
one to two orders of magnitude larger than the scalar coupling. The resonance frequencies 
and scalar couplings can be directly measured using multiple pulse sequences that average 



out the dipolar interaction, such as the MREV-8 sequence 19|, l20|, l2l|] . The intramolecular 
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dipolar coupling strengths for a partially ordered system are a result of the structure of the 



molecule as modified by the order parameters of the system 



18 



Quantum information 



processing has previously been explored in partially-oriented dipolar-coupled systems [23|, 
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25| . with as many as seven spins 26|], and strongly modulating radio frequency (rf) 
control sequences [27 1. 

We are interested in the control of two logical qubits that are encoded to protect against 
collective Oz noise. The logical subspace Sl for this encoding is the zero-quantum subspace 
of the Zeeman energy eigenstates 
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We describe the system in terms of the set of logical subspace Pauli operators 
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along with the nine bipartite terms such as cr^^a^^. Recently we reported on liquid state 
NMR experiments that demonstrate coherent control for a Bell state with this encoding 

we described a convenient subsystem pseudo-pure state [13] • Here we use a liquid crystal 
system to extend these studies, introducing a new symmetry for the spin system Hamilto- 
nian that leads us to expect that the logical encoding will be a more natural and efficient 
subspace for manipulating quantum information. In particular, the dipolar Hamiltonian has 



a portion that transforms as the exchange operator which has been shown to be convenient 
for subsystem encodings Q, 1^, Q, l3] ■ 

The goal of this work is to demonstrate three results: (1) improved quantum information 
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FIG. 1: 600 MHz proton spectrum of o-chloronitrobenzene (CNB) partially oriented by the liquid 
crystal solvent ZLI-1132, where the signal due to solvent protons has been suppressed by inserting 
a single cycle of the C-48 pulse sequence [29*] (r = 30/is) before the acquisition. The inset spectrum, 
collected stroboscopically under the MREV-8 sequence (r = 15^s), shows the four spins uncoupled 
with chemical shifts scaled by approximately 0.50. 

processing by using logical qubits, (2) an implementation of a DFS with dipolar-coupled 
spins, and (3) high fidelity control even when we have limited knowledge of the system 
Hamiltonian. The spin system used in these studies consists of the four protons of o- 
chloronitrobenzene (CNB) dissolved in Merck ZLI-1132 liquid crystal at 14 T field and a 
temperature of 300 K. The proton spins are strongly coupled to each other, and all of 
the resonances are not resolved in the 1-D NMR spectrum shown in Figure [H The order 
parameters for CNB aligned in a liquid crystal solvent have not been measured previously 
so we do not have this information on which to determine the internal Hamiltonian. We 
have made the following simple measurements to obtain reasonable estimates of the internal 
Hamiltonian: (1) a 1-D MREV-8 spectrum shown in Figured! and (2) 2-D correlation spectra 
between the chemical shifts under MREV-8 line narrowing and the full internal Hamiltonian 
(not shown). The MREV-8 spectrum indicates the chemical shifts, and the 2-D measurement 
provides a means of assigning the largest dipolar couplings to the appropriate chemical shifts. 
We measured the chemical shifts (in units of Hz) ui = 115, 1^2 = —234, z/3 = 204, and z/4 = 
—86 relative to an arbitrary transmitter frequency, and an incomplete set of approximate 
dipolar couplings (in units of Hz) di2 = —729, ^23 = —503, and ^34 = —1875. Although 
this limited description gives a very incomplete picture of the total system dynamics, it is 
sufficient for our purposes. 

We encoded two logical qubits into the four spin system where du provided the control 
elements to rotate the first logical qubit, ^34 similarly controlled the second logical qubit 
and (^23 controlled the interactions between qubits. To achieve the desired control fidelity we 
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found pulse sequences via the GRAPE (Gradient Ascent Pulse Engineering) algorithm [30 1. 
Starting with an initial guess of the control amplitudes for a pulse of specified length and 
number of intervals, the algorithm calculates the forward propagated states or unitaries at 
every time step, as well as the reverse propagated states or unitaries (from the desired target). 
It then uses the local gradient of the overlap between these forward and backward propagated 
states or unitaries, with respect to the control amplitudes to update the controls. Pulses 
were optimized for robustness to the unknown dipolar couplings (dis, du, ^24), uncertainties 
(±100 Hz) in the specified dipolar couplings {di2, (^23, (^34) and rf inhomogeneity. The three 
unknown dipolar couplings were set to a distribution of values (±100 Hz) centered about a 
"best guess" (in units of Hz) dis = 116, c/14 = —64, and ^24 = —170. All scalar couplings 
are small and were set to zero for pulse optimization. Although the most accurately known 
parameters in the internal Hamiltonian are the chemical shifts, control of the logical qubits 
did not rely on them 31|. 

The experimental goal was to create a pseudo-pure state over the logical qubits and then 
to entangle them in the form of a Bell state. We directly created the pseudo-pure state 
over the logical qubits via temporal averaging. This was accomplished in two steps. Under 
MREV-8 decoupling we prepared the states 
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relying on the differences in the chemical shifts of the four spins. To complete the pseudo- 
pure state preparation we used a GRAPE pulse to implement 
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We then implemented 
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which takes the input state vector |00)^ to the logical Bell state (|00)^ — lll)^) /V^- A 
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FIG. 2: The procedure and experimental results of creating a Bell state over two logical qubits 
encoded in the four dipolar-coupled protons of the CNB molecule. Beginning with the thermal 
equilibrium state, MREV-8 sequences along with a numerically optimized preparation pulse were 
used to create the pseudo-pure input state \00)j^ over the logical subspace. A numerically optimized 
entangling operation pulse converted the input state to the logical Bell state (|00)^ — |ll)^)/\/2- 
A set of 14 numerically optimized readout pulses were used to reconstruct the density matrices 
shown above. The experimentally measured logical input and Bell states have correlations of 0.90 
and 0.84 (respectively) with the numerically simulated states shown above. 

GRAPE pulse was found which performs this unitary operation over the logical subspace 
with a simulated fidelity of 0.99, accounting for coherent errors and the uncertainty in 
Hamiltonian parameters. In addition to the preparation and entangling pulses, we found 
fourteen readout pulses that transform every operator in the logical space into measureable 
operators, as in [32] . The full experiment is outlined in Fig. [21 

Information stored in the DFS respects a direct sum representation over the full Hilbert 
space. Ideal control operations over the logical qubits have the structure S = Sl^BSr, where 
Sl is the logical subspace of interest, and Sr is the remainder of the space. We attempted 
to prepare the initial state pj„ = |00)^ (00]^;^ © p^, where is arbitrary. The entangling 
operation over the full Hilbert space should also respect this symmetry Uent = Uent ® ^ent 
where we have complete flexibility in our choice of U^^. The challenge of course is that if 
the direct sum representation is not maintained, the logical information can be corrupted 
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or can leak out of the logical subspace [28|, l33|, [SJ] . 



Uprep and Uent Were optimized over the range of dipolar couphngs and rf inhomogeneity 
described above. The goodness function maximized by the GRAPE algorithm was the logical 
subspace fidelity [I4], which includes a penalty for pulses that permit leakage. was 
chosen arbitrarily. Finally, the readout pulses necessarily operated over the entire Hilbert 
space. They were simply designed to efficiently transform selected logical operators into 
observables and to be robust over the dispersion of coupling constants and rf inhomogeneity. 

The reconstructed density matrices over the logical degrees of freedom are shown in Fig. 
O As expected, most of the observed errors arise from the initial state preparation and the 



readout sequences. The normalized state correlations over the logical subspace 



]J are: 



Corr {ptr^pTJ) =0.90 (10) 
Corr(p-,p^^S,)=0.84 (11) 

The normalized correlations between the experimentally measured states and the ideal states 
are of course lower since they include more of the errors due to state preparation and readout 

Con {pfr',pZ^) =0.83 (12) 
Corr(p^^:;]',pa=0.76 (13) 

We tested the robustness of our control over the logical subspace by calculating the 
average state correlation that would be expected as the dipolar frequencies are varied. The 
average correlation between simulation and experiment was found to decrease very little (less 
than 0.03) even when all six dipolar frequencies were simultaneously varied by ±100 Hz. 
This is expected, since the pulse sequences were engineered to be robust to these variations 



351]. This result shows that even in the presence of some uncertainty in our Hamiltonian 
parameters, we have reliable control over the logical subspace of our four proton system. 
Finally we compare our control over the encoded logical qubits to that we have over the 
individual spins. In numerical simulations, we compared the logical qubit entangling pulse 
implemented in the experiment to the best spin qubit entangling pulses found by the same 
pulse optimization methods, optimizing over the full Hilbert space. The decoherence rate 
of the logical qubits is slower than the decoherence rate of the spin qubits due to the noise 
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protection of the DFS. The spin qubit control sequences must be long enough to resolve 
the differences in spin resonance frequencies (|z/j — Ukl) and therefore are necessarily longer 
in time than the logical qubit control sequences. As a result, when the output states are 
attenuated by decoherence, the fidelity of entangling logical qubits is approximately 0.90, 
while the fidelity of entangling spin qubits is significantly lower, approximately 0.70 in the 
best case js^ . This result quantifies the improvement in control enabled by our logical qubit 
encoding. 

In conclusion, we have shown that improved coherent control is achieved by encoding 
logical qubits in our system of four dipolar-coupled protons. Control in this system was 
demonstrated experimentally by creating a pseudo-pure state in a DFS and applying a 
unitary transformation to create a logical Bell state. Though the inital and final states in 
our experiment are not pure, the implemented transformation is still entangling. Some of 
the system Hamiltonian parameters are imprecisely known, and high fidelity control was 
achieved by engineering pulse sequences that are robust to these uncertainties. The dipolar 
couplings in the system Hamiltonian, in the strong coupling regime, provides a natural 
setting for experimental studies of logical qubit encodings 23|]. In the future, liquid crystal 
solvent NMR QIPs could be used to explore more complex logical encodings in larger Hilbert 
spaces. In these larger systems we expect that the nearest-neighbor couplings will not 
dominate to such an extent, and that the system will not map onto a one- dimensional chain, 
as did this 4-spin example. 
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ported under an NSF Graduate Research Fellowship, in addition to support by the National 
Security Agency (NSA) under Army Research Office (ARO) contract W911NF05 10469 and 
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